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Abstract 

An overview of NASA's ongoing efforts to 
develop an airfoil icing analysis capability is 
discussed. An indication is given to the 
approaches being followed to calculate ( 1 ) the 
water droplet trajectories past the airfoil, ( 2 ) 
the buildup of Ice on the airfoil, and (3) the 
resultant changes in aerodynamic performance due 
to the leading edge ice accretion. Examples are 
given of current code capabi J ? ties/ 1 imi tat ions 
through comparisons of predictions with experimen- 
tal data gathered in various calibration/validation 
experiments. A brief discussion of future efforts 
to extend the analysis to handle three-dimensional 
components is included. 

Li st of Symbol s 
c airfoil chord 

C(j airfoil drag coefficient 

C) airfoil lift coefficient 

Cmc /4 airfoil moment coefficient 
C p airfoil surface pressure coefficient 

LWC icing cloud liquid water content 

MVD Icing cloud volume median droplet diameter 
S surface distance 

Ta, freestream total temperature 

U local velocity 

Ue edge velocity 

Vqo freestream velocity 

x, y coordinates 

x^ r transition location 

a airfoil angle-of-attack 

fi local water droplet collection efficiency 

9 angular position 

Introduction 

One of the major elements of the NASA Air- 
craft Icing Research Program is the development 
and validation of a series of mutually compatible 
computer codes which predict the details of an 
aircraft Icing encounter. These codes once vali- 
dated can be used for: 

(1) Making preliminary design studies to 
ascertain component sensitivity to icing. 


(2) Making performance predictions of pro- 
posed ice protection systems. 

(3) Conducting computer based certif ication/ 
qualification studies to reduce the amount of 
required icing flight testing. 

(4) Providing more accurate and realistic 
icing effects inputs for use in flight training 
simulators . 

Figure lisa flowchart which shows the many 
codes required to form such an icing analysis 
methodology as well as the codes currently under 
development by NASA. Reference 1 provides an 
overview of the icing analysis activities (analyt- 
ical and experimental) underway. 

The purpose of this paper is to review the 
progress on one element of the overall activity, 
namely the unprotected airfoil icing problem. 

Simply stated, the problem being addressed is to 
develop/validate computer codes which predict ( 1 ) 
the buildup of ice on an unprotected airfoil and 
( 2 ) the resultant aerodynamic degradation levels 
due to that ice accretion shape. A validated air- 
foil icing analysis capability Is felt to be the 
core capability for the overall icing analysis 
methodology shown in Fig. 1. This airfoil icing 
analysis capability will be extended to handle 
three-dimensional components (e.g., swept wings) 
and eventually complete aircraft configurations. 

Figures 2 and 3 pictorially show the problem 
being addressed. A wing section of constant chord 
was placed in the NASA Lewis Icing Research Tunnel 
(IRT) and subjected to specified icing conditions 
(cloud liquid water content, average droplet diam- 
eter, free stream temperature and icing time) and 
aerodynamic conditions (velocity and angle-of- 
attack). The figure shows the resultant leading 
edge ice accretion in overall and closeup views. 
The rough Irregular nature of the ice accretion 
should be noted- 

Figure 3 shows in a generalized fashion the 
decrease in aerodynamic performance which can 
occur due to leading edge ice accretions; namely: 

(1) An Increase in wing section drag even at 
low angles-of-attack 

(2) A decambering of the airfoil due to a 
thickened upper surface boundary layer 

(3) A reduction in Ci max and, premature 
stall due to separation of the airfoil upper sur- 
face boundary layer. 

For any given airfoil geometry, the ice 
accretion can change dramatically in cross sec- 
tional shape as the environmental and aerodynamic 
parameters are varied. The level of severity of 
aerodynamic performance degradation will also vary 
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significantly depending on the ice accretion cha- 
racteristics as well as the basic airfoil perform- 
ance characteristics. 

Previous experimental studies have determined 
many of the key aspects of the airfoil icing prob- 
lem. These aspects must be adequately modeled/ 
predicted in any airfoil icing analysis methodology 
(Fig. 4). While an airfoil ice accretion shape is 
typically small relative to the airfoil size (usu- 
ally only a few percent chord in maximum dimension), 
its effect on the airfoil flowfield and thus the 
airfoil aerodynamic performance can be profound. 

The ice accretion shape is highly rough and irreg- 
ular (as the closeup portion of Fig. 2 suggests). 

The roughness elements (size and distribution) are 
much different than the classical sand grain ele- 
ments used in wind tunnel test programs over the 
years. The developing boundary layers on both the 
suction and pressure sides will tend to separate 
in the near vicinity of the ice accretion. At 
lower angles-of-attack, boundary layer reattach- 
ment from either or both surfaces will usually 
occur aft of the ice accretion, but the reattached 
viscous layer(s) will be thickened and distorted. 

As the angl e-of-attack is changed, the relative 
sizes of the suction and pressure surface separation- 
reattachment zones will change. In particular as 
the angle-of-attack is increased, the suction sur- 
face zone will grow in size until no reattachment 
occurs and airfoil stall results. Airfoil stall 
can occur several degrees below the stall angle 
observed for clean airfoils. The increase in dray 
occurring at the lower angles-of-attack is due to 
thickened, distorted boundary layers which develop 
over the airfoil aft of the leading edge ice 
accretion. 

The mutually compatible computer models 
required to form an airfoil icing analysis capa- 
bility are shown in Fig. 5. The overall approach 
to the problem will first be discussed and then 
each module will be discussed. 

Overall Approach 

As Fig. 5 Indicates, a potential flow analy- 
sis is first performed to determine the inviscid 
flowfield around the clean airfoil. The trajectory 
module computes the paths of water droplets past 
the airfoil and determines the rate at which drop- 
lets are impacting the airfoil surface at each 
point. This represents the local rate of water 
deposition on the airfoil. 

The ice accretion module calculates the local 
rate of ice growth around the airfoil. The dis- 
tribution of ice growth rates around the airfoil 
surface allows an Ice shape to be predicted once a 
user specified time Increment (called the time 
step) Is selected. This time Increment is chosen 
such that the effect of ice buildup on airfoil 
flowfield and thus droplet trajectory paths is 
thought to be minimal. 

Once the ice shape is calculated, then a vis- 
cous flowfield calculation is performed using the 
predicted ice shape If aerodynamic performance 
degradation levels are required. 

If It is desired to continue to computation- 
ally accrete Ice on the airfoil, then the looping 
process is repeated for as many time increments 


as is required to reach the overall icing encoun- 
ter time. Each successive loop begins with a 
potential flowfield analysis of the airfoil with 
the leading edge modified by the ice accretion 
computed from the previous step(s). 

It is highly desirable for the separate 
inviscid and viscous flow analyses to be combined 
(i.e., replace the inviscid analysis step with 
the viscous analysis). However, a significant 
increase in computational resources currently is 
required to complete the viscous analysis (minutes 
on a Cray XMP) versus that required to complete 
the inviscid analysis (minutes on an IBM 370). 
Obviously as the ice shape grows in size and domi- 
nates the airfoil leading edge flowfield, viscous 
effects (boundary layer separation/reattachment) 
will become so dominant that the simplified invis- 
cid analysis will no longer be appropriate. 

The following sections will look at the 
modules in more detai 1 . 

Inviscid Flowfield/Drpplet Trajectories 

The flowfield about the airfoil (dean and 
iced) is calculated using a potential flow 
analysis - a second order panel code developed by 
Hess-Smith. Flowfield information - individual 
panel geometry and source/vor t i c i ty strengths - is 
passed to the trajectory module to allow calcula- 
tion of the water droplet paths past the airfoil. 
The droplet equations of motion are integrated 
using a predictor-corrector scheme optimized for 
stiff systems of equations. Flowfield velocity 
values at various points on the airfoil surface 
required for integration of the equations of 
motion are determined by summing the individual 
contributions of all panels. 

The accuracy of the trajectory analysis is 
being assessed by comparison with an experimental 
droplet impingement data base being developed 
jointly by NASA and the FAA.2 Selected compari- 
sons are shown in Fig. 6 for data gathered on a 
clean cylinder model as well as on an "iced" cyl- 
inder model. The figure also shows droplet col- 
lection efficiency predictions made using a 
Navi er-Stokes flowfield analysis in place of the 
Hess-Smith potential code. Both sets of predic- 
tions generally agree well with the clean cylinder 
data although the tails of the collection effi- 
ciency curves are broader for the data than for 
the predictions. This was expected as the calcu- 
lations were made for a single droplet size (the 
average diameter) while the actual droplet cloud 
had a spectrum of sizes. 

The agreement is not as good for the "iced 11 
cylinder, especially in the so-called horn 
regions. Experimental measurements of the surface 
pressure distribution about the ice shape indicated 
the boundary layer separated just aft of each 
horn. Since potential flow analyses do not 
account for separation, the local velocities in 
the vicinity of each horn are predicted to be 
larger than actually exist. This inadequacy in 
the flowfield modeling was felt to be reason for 
the overprediction of the collection efficiency 
levels in the vicinity of the horns and led to 
the trajectory calculations using the Navier- 
Stokes analysis. Suprisingly the inclusion of the 
viscous effects did not significantly decrease 
the collection efficiency levels near the horn. 
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The reason for this discrepancy is not known but 
it may suggest some inaccuracy in the experimen- 
tal approach. The surprising level of agreement 
in collection efficiency predictions using the 
panel and Navi er-Stokes codes for both the clean 
and iced cylinders suggests the panel code calcu- 
lations can provide reasonable estimates of the 
collection efficiency distributions for clean and 
iced ai rfoi 1 s . 

As already stated, in the long term it is 
desirable to replace the inviscid analysis with a 
viscous analysis in the airfoil icing analysis 
methodology to properly account for boundary layer 
separation - reattachment. 

Ice Accretion 

In spite of the fact that aircraft icing has 
been studied since the 1920's, the understanding 
of the governing physical processes of ice accre- 
tion is limited. Consequently, ice shapes are 
predicted using a control volume approach as 
depicted in Fig. 7. Mass and energy balance equa- 
tions are formulated accounting for the effects 
shown in the figure. The equations are solved to 
determine the fraction of incoming water which 
freezes In each control volume. Any water which 
does not freeze in a control volume is assumed to 
flow back to the immediately aft control volume 
(so called runback water). Required inputs from 
the flowf ield/trajectory analysis include the 
collection efficiency and surface velocity 
distributions . 

Figure 8 shows representative comparisons of 
predicted versus measured ice shapes. In each 
case, the ice shape prediction is compared to an 
ice shape grown in the NASA IRT. The colder tem- 
perature condition resulted in a rime ice condi- 
tion while the warmer temperature resulted in a 
glaze ice accretion. (Glaze ice shapes tend to 
have horn characteristics which result In large 
levels of aerodynamic degradation). Also note 
that the figure shows the predicted ice shape(s) 
were grown in time steps, each of 30 seconds dura- 
tion. The agreement in predicted versus measured 
ice shape for both cases is judged to be accepta- 
ble. Currently efforts are underway to make a 
large number of Ice shape comparisons with not 
only data gathered in many icing ground tests 
facilities but also with natural icing flight 
data. This extensive set of comparisons will 
serve as a code calibration/validation effort for 
the ice accretion analysis. 

A significant input to the ice accretion cal- 
culation procedure has been shown to be the con- 
vective heat transfer coefficient (h) distribution. 
Currently, the h distribution is calculated using 
an Integral boundary layer approach described by 
White-* and Kays and Crawford. 4 An experimental 
heat transfer coefficient data base has been 
acquired 5 to evaluate the adequacy of the predic- 
tion approach, and a comparison Is shown in 
Fig. 9. The cylinder Ice shape shown is the same 
as the shape used in the trajectory impingement 
studies. The major disagreement is seen to be In 
the vicinity of the horn and downstream of the 
horn. The disagreement is attributed to the 
potential flow prediction of surface velocity dis- 
tribution not accounting for the observed boundary 
layer separation. 


As already Indicated (Fig. 2), Ice accretion 
shapes can exhibit a large degree of surface rough- 
ness and an attempt has been made to account for 
this effect in the ice accretion analysis. A dis- 
tribution of sand grain roughness is input and the 
heat transfer rates are increased according to the 
predictions of the integral boundary layer approach. 
The heat transfer data base acquired included 
testing ice shapes to which sand grains of known 
size were affixed. One prediction/experimental 
comparison of this rough surface heat transfer is 
shown in Fig. 10. The same general remarks can be 
made as were already made for the smooth surface 
ice shape. Also a comparison of Figs. 9 and 10 
reveals that the increase in measured heat trans- 
fer due to roughness was predicted reasonably well 
by the integral boundary layer routine. 

Experimental efforts are currently being con- 
ducted to acquire additional heat transfer data 
for both smooth and rough surface airfoils for 
comparison with predicted levels. The data base 
includes measurements of the effect of free stream 
turbulence on heat transfer levels. Also, efforts 
are underway to determine more accurate ways of 
characteri zi ng the surface roughness than the sim- 
ple sand grain element approach currently used. 

In order to improve the physical understand- 
ing of the ice accretion process, close-up photo- 
graphs and movies have been taken of the icing 
process on an airfoil. 5 Such data has been 
acquired both in artificial and natural icing con- 
ditions. These photographic results suggest that 
some of the assumptions contained in the control 
volume model of ice accretion may be in error. 

In particular, the model allows for continuous 
flow of water on the surface. That is, any water 
not frozen in a control volume is allowed to run 
back to the adjacent control volume. The photo- 
graphic evidence suggests that after an initial 
transient, surface water (if it exists) does not 
appear to move. A multi zone model that has been 
proposed to account for this discrepancy is 
depicted in Fig. 11. Surface water is allowed to 
run in Zone A which would be characterized as a 
smooth region. Zone C would be a transition 
region which would have a varying roughness while 
Zone B would be a region which would have no run- 
back. The transition region (Zone C) would move 
toward the stagnation point as the icing time 
Increases and thus Zone A would decrease in size 
with icing time. Currently, the position of Zone 
C is determined empirically but it is most desir- 
able to be able to relate it to the governing 
physics. Currently this multizone model is being 
added to the ice accretion module to evaluate the 
resultant improvement in Ice shape predictions. 

Also, efforts are underway to model some of 
the other aspects of the ice accretion process 
noted in the photographic data that are not 
included in current ice accretion analysis. These 
effects include droplet impact and possible splash- 
ing and water beading into large droplets. Addi- 
tional efforts (experimental and analytical) are 
underway to better understand the water film 
dynamics effects for inclusion into the ice accre- 
tion modeling. It is hoped that inclusion of 
these physical effects in the ice accretion 
analysis will contribute to improved ice shape 
predictions. 
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Aerodynamic Performance 

Two analysis codes are being applied to pre- 
dict airfoil aerodynamic performance degradation 
due to ice accretion: ( 1 ) a Reynolds averaged 

Navier Stokes code' (ARC2D) and (2) an interactive 
boundary layer code (IBL ). 8 Both of these codes 
were originally developed to handle "clean" geome- 
tries and are being extended to handle iced air- 
foils. In order to properly evaluate these codes, 
a comprehensive validation data base is required. 
The initial approach being followed to acquire 
such a data base Is summarized in Fig. 12. A 
NACA001 2 airfoil model was modified to have a lead- 
ing edge with an idealized ice shape. This ice 
shape was chosen to have the gross cross sectional 
features of an ice shape grown in the IRT but to 
have a well defined surface description to allow 
inputting into the flow analysis codes. Detailed 
surface pressure distributions and wake survey 
probe measurements allowed global force and moment 
performance to be documented. Boundary layer 
velocity profiles on both suction and pressure 
surfaces were acquired using both pitot tubes and 
single element hot film probes. The profile mea- 
surements were concentrated in the separation- 
reattachment zones. Flow visualization data was 
acquired to complement the quantitative data base. 
Reference 9 provides more details of the experi- 
mental program. 

Figure 13 summarizes the predictions of the 
ARC2D and IBL codes when compared to the experi- 
mental data base already discussed. The predic- 
tions of both codes for lift and drag coefficient 
levels for the lower angles-of-attack agree well 
with the data. At the higher angles-of-attack, 
the IBL code tended to underpredict the measured 
drag levels. At these higher angles-of-attack, 
the ARC2D code predicted that an unsteady flow- 
field existed (i.e., no steady solution was 
reached when the code was run in the time accu- 
rate mode). The ARC 2 D results for those higher 
angles-of-attack were determined by averaging the 
surface pressure distributions over one cycle and 
then calculating the resultant C] , Cmc/ 4 , and 
Cq values. The IBL code also had some difficulty 
achieving a converged solution at these conditions 
but due to the formulation of the interactive 
boundary layer approach, no such averaging proce- 
dure was possible. 

While the global comparisons of code predic- 
tions of Ci and Cq versus experimental mea- 
surements yielded generally favorable results, 
significant differences did occur when comparisons 
were made of surface pressure distributions and 
velocity profiles. Both codes tend to predict 
upper surface boundary-layer reattachment to occur 
prior to the measured reattachment points. 

Figure 14 shows some selected velocity profile 
comparisons. Note that both codes predicted sepa- 
rated profiles at X/C * 0.005 which were signifi- 
cantly different than the measured profile . The 
ARC2D code predicted reattachment to have occurred 
at station 0.043 but the experimental profile sug- 
gested that reattachment had not occurred. Also, 
the predicted static pressure distributions showed 
much greater negative pressure peaks near separa- 
tion than were observed to occur in the experiment. 


Three possible sources of these discrepancies 
were noted and are currently being investigated: 
( 1 ) ice shape def i n 1 t ion/f lowf i el d grid genera- 
tion, ( 2 ) boundary layer transition, and ( 3 ) tur- 
bulence model ling. 

Figure 15 shows the results of a grid defini- 
tion study conducted with the ARC 2 D code using an 
ice accretion shape predicted using the ice accre- 
tion analysis module already discussed. Note the 
Ice shape is not smooth and well defined like the 
idealized ice shape already discussed. Two NASA 
Ames developed grid generation codes were used in 
the study: a Poison solution approach (GRAPE) and 

a hyperbolic generation approach (HYPGRID) . The 
curvilinear body fitted grids were 253 (along the 
surface) X64 (into the flowfield). 

The GRAPE code was used with points uniformly 
distributed around the surface while the HYPGRID 
code was used with both uniform distributions and 
densely packed leading edge distributions. Also, 
the HYPGRID code was run with more points placed 
closer to the airfoil surface so as to capture 
more details of the growth of the viscous layers. 

A comparison of the ARC2D predictions on a C] , 

Cq basis indicated the only appreciable differ- 
ences noticed were in the C] predictions for the 
higher angles-of-attack. Packing points in near 
the ice surface (both along the surface and nor- 
mal to the surface) resulted in delayed stall pre- 
diction and thus higher predicted C] . This is 
attributed to improved computations of vorticity 
generation and resolution of the boundary layer 
growth . 

The results of this initial study failed to 
show any significant differences when using grids 
from the GRAPE and HYPGRID codes. These studies 
are continuing. Also ice shape sensitivity stud- 
ies are being conducted with the IBL code. Since 
the IBL analysis combines potential flow and bound 
ary layer routines, no grid generation is required 
These studies are aimed at determining how 
detailed the input ice shape must be to yield 
accurate results from the aerodynamic analysis. 

Previous clean airfoil studies 10 have identi- 
fied that boundary layer transition location can 
have a significant effect on airfoil performance 
predictions. A limited study was conducted with 
the ARC2D code to assess the importance of transi- 
tion location specification for the iced airfoil. 
Figure 16 shows the results for the computations 
of the flow about the idealized ice shape previ- 
ously discussed for two specific transition loca- 
tions. The forward location corresponded to the 
tip of the upper horn while the aft location cor- 
responded to the point midway between the separa- 
tion and reattachment locations. There was no 
effect on the computed values for lift and drag 
but a noticeable effect did occur for the velocity 
profiles in the separated flow zone and the size 
of the separated flow bubble. 

The aft specification of the transition loca- 
tion resulted in a larger bubble with a reattach- 
ment point which agreed more closely with the 
experimental measurements. However, the computa- 
tions forced a rel ami ni zat ion as the fluid moved 
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upstream forward of the transition point. Conse- 
quently, the magnitude of the reverse flow velocity 
is overpredicted. For the forward specification 
of transition, the separated flow zone was com- 
pletely turbulent and the velocity levels near the 
surface were in better agreement with the experi- 
mental data. 

The ARC2D calculations shown in the previous 
figures were accomplished using a Baldwi n-tomax 
turbulence model. Recently, investigations have 
been initiated to look at the effect of alternate 
turbulence models on iced airfoil flowfield pre- 
dictions. The initial results of that study are 
summarized in Fig. 17. Two turbulence models have 
been used for comparison with the Baldwin-Lomax 
results: the Johnson-King model (an ordinary dif- 

ferential equation model) and the k-c model (a 
two equation model). A comparison of the leading 
edge surface pressure distributions shows that no 
noticeable effect could be seen. All three sets 
of calculations resulted in a significant overpre- 
diction of the suction peak and no prediction of 
the observed pressure plateau. In addition, all 
three calculations indicated a smaller separation 
zone than was measured. The comparison of the 
predicted velocity profiles in the separation zone 
as well as the global values of C] , C<j, and Cmc /4 
also failed to show any significant differences. 

These results suggest that implementation of 
a more sophisticated turbulence model in the iced 
airfoil analyses is not warranted at this time. 

The additional CPU time required to perform the 
calculations does not appear to yield any improve- 
ment In results over those acquired using a simple 
algebraic equation turbulence model. Efforts are 
continuing to look at other algebraic equation 
models which are available. These results also 
suggest that a turbulence model “tuned" to handle 
separated-reattached flows will most likely be 
required . 

Concluding Remarks 

This paper has reviewed the current status 
of an ongoing NASA activity to develop and vali- 
date an airfoil icing analysis capability . The 
selected results shown indicate that first genera- 
tion analysis codes exist to predict ( 1 ) ice 
accretion shapes and ( 2 ) resultant aerodynamic 
performance levels. The codes have not as yet 
been integrated Into a single analysis methodology 
due to the large computer times required to per- 
form the aerodynamic analysis, but this is a 
desired longer term goal. Work is underway to 
improve the understanding of the key physics of 
the ice accretion process so that the ice accre- 
tion analysis code can be modified to improve the 
prediction of ice shapes. Two aerodynamic analy- 
sis codes are being evaluated and both show prom- 
ise for predicting airfoil degradation due to ice 
accretions. Limitations in the two aerodynamic 
analysis have been identified to include: ice 

shape description/grid generation, transition 
specification, and turbulence modelling. Efforts 
are underway to address these deficiencies. 

An additional limitation which also must be 
added to both the ice accretion and aerodynamic 
analyses is an adequate treatment of the extreme 
levels of surface roughness which can exist on ice 


shapes, particularly those formed at warmer tem- 
peratures (close to freezing). As already indi- 
cated these roughness elements are much different 
than the classical sand grain roughness elements 
used in wind tunnel testing. It is felt that 
these roughness elements must be modelled since 
they can have a significant effect on the boundary 
layer growth and thus on heat transfer rates and 
aerodynamic degradation. 

A validated airfoil Icing analysis methodol- 
ogy will provide a most useful computational tool 
to the icing community which can be used for a 
number of purposes. The methodology will also 
provide the basis for developing a three- 
dimensional icing analysis methodology to handle 
geometries such as swept wing configurations. 

Figure 18 suggests the distinct three-dimensional 
character that ice accretion shapes can acquire 
for near freezing temperature conditions. Cur- 
rently efforts are underway to evaluate various 
three dimensional Navier-Stokes codes for possi- 
ble application to the three-dimensional icing 
problem. Also, a three-dimensional version of the 
interactive boundary layer code is under develop- 
ment and will be employed. An experimental pro- 
gram has been initiated to measure degradation of 
wing aerodynamic performance due to a known ice 
shape. Tests will be conducted for both straight 
and swept wing conf i gurations . Detailed flowfield 
measurements will be made using techniques similar 
to those used during the airfoil icing experiments. 

Clearly, the analysis of airfoil icing is a 
difficult problem for the three-dimens ional meth- 
ods. The three-dimens ional problem will provide 
a much greater challenge to the icing analysis 
community. The development of such a capability 
will allow evaluation of civil and military air- 
craft configurations performance in icing condi- 
tions. Such an analysis capability is considered 
to be valuable to the aerospace community and is 
a desired end goal of the NASA aircraft icing 
analysis program. 
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FIGURE 1. - AIRCRAFT ICING ANALYSIS METHODOLOGY. 
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FIGURE 2. - THE AIRFOIL ICING PROBLEM. 




DECAMBER I NG 

FIGURE 3. - AERODYNAMIC PERFORMANCE DEGRADATION DUE TO ICING. 
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FIGURE A. - KEY ASPECTS OF AIRFOIL ICING, 



FIGURE 5. - STEPS IN AIRFOIL ICING ANALYSIS METHODOLOGY. 
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FIGURE 8. - COMPARISON OF ICE SHAPE PREDICTIONS WITH AIRFOIL ICING DATA. 



FIGURE 9. - COMPARISON OF CALCULATED AND MEASURED HEAT 
TRANSFER COEFFICIENTS ON SMOOTH SURFACE "ICED" 
CYLINDER MODEL. 



CYLINDER ANGLE, 0 

FIGURE 10. - COMPARISON OF CALCULATED AND MEASURED HEAT 
TRANSFER COEFFICIENTS ON ROUGH SURFACE "ICED" CYLINDER 
MODEL. 
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FIGURE 11. - MULTIZONE MODEL OF ICE ACCRETION PROCESS. 





FORCE AND MOMENT DATA 
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FIGURE 12. - CODE VALIDATION DATA BASE FOR ICED AIRFOIL PERFORMANCE. 




FIGURE 13. - COMPARISON OF ICED AIRFOIL CODE PREDICTIONS WITH EXPERIMENTAL 
MEASUREMENTS. 
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GRID SPACING 

ALONG NORMAL 

SURFACE TO SURFACE 

O DENSE L.E. NORMAL 

O EVEN NORMAL 

v DENSE L.E. DENSE 

a EVEN DENSE , 

D EVEN NORMAL 


GRID CODE 


HYPGRID 


GRAPE 



FIGURE 15. - RESULTS OF GRID GENERATION STUDY. 
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X/C 



U/U e 


x TR /c = 0.05 C£ = 0 c d = 0.302 

FIGURE IS. - EFFECT OF BOUNDARY LAYER TRANSITION SPECIFICATION IN NAVIER-STOKES 
PREDICTED VELOCITY PROFILES IN SEPARATION-REATTACHMENT ZONE, Of = 0. 



0 .125 .250 0 .125 .250 0 .125 .250 

x/c x/c x/c 


BALDWIN-LOMAX JOHNSON-KING K-E MODEL 

MODEL MODEL 

FIGURE 17. - EFFECT OF TURBULENCE MODEL ON NAVIEP-STOKES PREDICTIONS OF LEADING EDGE PRESSURE 
DISTRIBUTIONS, a = 0°. 
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FIGURE 18. - THE EFFECT OF WING SWEEP ON ICE ACCRETION 
CHARACTERISTICS. 
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